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INTRODUCTION 


Two  of  the  primary  pieces  of  the  diagnostic  information  in  pulsed-power 
experiments  are  the  voltage  and  the  current  provided  by  the  pulser  at  the 
load.  In  many  instances,  particularly  with  voltages  on  vacuum  devices,  it  is 
not  possible  to  make  measurements  as  close  to  the  load  as  one  would  like. 
Often  an  inductive  correction  is  used  to  relate  the  voltage  sensed  at  one 
position  to  the  actual  voltage  at  the  position  of  interest.  However,  such  a 
correction  is  only  an  approximation  with  limited  validity. 

A  better  method  to  infer  the  voltage  and  current  at  some  position  from 
data  measured  at  another  position  uses  a  transmission-line  circuit  and  a 
transmission  matrix  approach.  There  are  three  requirements  for  the  applica¬ 
tion  of  this  method:  first,  the  voltage  and  current  must  be  sensed  at  a 
common  position;  second,  the  record  of  these  data  must  include  information 
before  and  after  the  time  frame  of  interest  by  a  time  equal  to  the  transit 
time  from  the  sensors  to  the  calculation  position;  and,  third,  the  geometry 
connecting  the  calculation  position  with  the  data  position  must  be  character¬ 
ized  by  a  system  of  transmission  lines. 

One  necessity  for  the  practical  implementation  of  this  method  is  a  comput¬ 
erized  system  for  pulsed  digitization  and  data  processing.  With  the  more 
widespread  use  of  high-speed  transient  digitizers,  this  otherwise  awkward  (and 
therefore  little-used)  method  is  quite  easy  to  implement. 


2.  ANALYSIS 

For  our  purposes  we  consider  that  the  region  connecting  the  field  sensors 
with  the  position  of  interest  can  be  described  by  a  sequence  of  lossless 
transmission  lines  each  described  by  an  impedance,  Z^,  and  a  one-way  transit 
time,  T^. 

In  each  transmission  line  there  can  be  a  wave  propagating  to  the  right  and 
to  the  left,  and  each  wave  has  a  frequency -independent  speed.  For  the  right¬ 
going  wave  of  voltage,  V,  the  current  at  a  given  point  in  the  line  is  given 
by  V/Z,  while  for  the  left-going  wave,  V,  the  current  is  -V/Z.  Consequently, 
the  total  curjent^  I,  is  given  by  I  =  (v  -  v)/Z,  whereas  the  total  voltage  is 
given  by  V  =  V  +  V.  Combining  the  expressions,  one  can  express  the  traveling 
waves  in  terms  of  I  and  V; 


V 


-  2(V 


+  ZI) 


and 


(1  ) 


V  =  I(V  -  ZI)  . 

We  assume  we  know  as  functions  of  time  the  voltage,  V^,  and  current,  1^, 
to  the  left  of  the  ith  transmission  line  (fig.  1  ).  What  we  want  to  determine 
are  the  functions,  V^+1  and  I^+1 ,  which  describe  the  voltage  and  current  to 
the  right  of  the  transmission  line. 
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Figure  1.  Definition  of  the  parameters:  and  I1  are  the  total 

voltage  and  current  to  the  left  of  the  ith  transmission  line;  and 

I^+1  are  the  total  voltage  and  current  to  the  right;  and  are, 
respectively,  the  right-  and  left-going  voltage  waves  in  the  ith 
transmission  line  at  its  left  end;  and  $R  and  trR  are,  respectively,  the 
right-  and  left-going  voltage  waves  in  the  same  line  at  its  right  end. 


At  a  given  time,  t,  the  voltage  and  current  to  the  right  of  the  transmis¬ 
sion  line  are  expressed  by 

vi+1 (t)  =  VR(t)  +  VR(t) 


and 


(2) 


ziIi+1 (t)  =  VR(t)  -  VR(t)  , 


in  which  the  subscript  R  identifies  the  traveling  waves  as  they  appear  in  the 

right  end  of  the  ith  transmission  _line.  ^These  values,  V  and  VR,  at  time  t 

are  equal  to  the  left-end  values,  V  and  V  ,  at  times  offset  by  the  transit 

.  tl  Li 

time,  : 

Vt}  =  , 

<-  «-  (3) 

VR(t)  “  * 

Using  equation  (1),  we  can  express  the  time-offset  traveling  voltages  at  the 
left  side  in  terms  of  the  total  voltage  and  current  measured  there: 

VL(f-Ti)  ={vi(t-Ti) 

and  ( 4 ) 

VL(t+Ti)  =^Vi(t+Ti)  -  J  ZiIi(t+Ti)  . 


We  now  combine  equations  (4),  (3),  and  (2)  and  obtain  the  expressions  desired 

for  and  Ii+<j(t)  in  terms  of  the  known  functions,  and  1^: 
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Vi+1<t)  =^[vi(t-Ti)  +  ZiIi(t-Ti)  +  Vi(t+Ti)  -  ZiIi(t+Ti)j  , 
Ii+1(t)  =l[vi(t-Ti)  +  ZiIi(t-Ti)  -  Vi(t+Ti)  +  ZiIi(t+Ti)]/Zi 


(5a) 


(5b) 


Equations  (5a)  and  (5b)  embody  the  result  of  the  transmission  operation 
through  one  element  of  the  circuit  and  are  the  equations  used  in  the  computer 
application  of  this  method. 

For  geometries  that  can  be  described  by  a  series  of  transmission  lines  of 
a  variety  of  impedances,  equations  (5a)  and  (5b)  can  be  used  to  step  from  one 
junction  to  the  next  until  the  voltage  and  current  are  derived  for  the  posi¬ 
tion  of  interest.  In  carrying  out  this  procedure  neither  V  nor  I  nor  the 
geometry  need  be  specified  to  the  left  of  the  point  at  which  V  and  I  are 
measured  or  to  the  right  of  the  position  of  interest.  Consequently,  derived 
data  can  be  used  directly  to  indicate  the  load  characteristics. 

According  to  equations  (5a)  and  (5b),  to  derive  or  Ii+1  at  a  time,  t, 

the  given  and  1^  must  both  be  known  at  times  t  +  and  t  -  T^ .  This 

requirement  chains  itself  across  a  succession  of  transmission  lines  so  that 
the  derived  data  are  less  extensive  than  the  measured  data  at  each  end  of  the 
trace  by  the  sum  of  the  individual  T^. 

3.  RELATION  TO  OTHER  DATA -CORRECT I ON  FORMULAS 

The  result  contained  in  equations  (5a)  and  (5b)  should,  with  the  proper 
assumptions,  simplify  giving  the  inductive  correction  term  LI  if  certain 
simplifying  assumptions  can  be  made.  To  investigate  what  these  assumptions 
might  be,  we  write  equation  (5a),  grouping  together  the  voltage  and  current 
terms : 

Vi+l(t)  =l[vi(t-Ti)  +  Vi(t+Ti)J  -  TiZi[li(t+Ti)  -  ( t-Ti)  ]  /( 2Tj_)  .  (6) 

The  factor  T^Z^  seen  above  is  identically  the  inductance  L^  of  the  ith  trans¬ 
mission  line.  If  the  variation  in  time  of  is  approximately  linear  from  t  - 
T^  to  t  +  Ti,  then  1  /  2[  Vf(  t-TjJ  +  V^(t+Tj1]  approximates  V^(t).  Furthermore, 
if  the # same  condition  holds  for  1^,  then  [l^(t+T^)  -  I^(  t-T^J  ] /( 2T^)  approxi¬ 
mates  I^(t).  If  both  and  1^  are  approximately  linear  over  t  -  T-  to  t  + 
T^,  we  see  that  equation  (6)  reduces  to 

vi+i(t)  -  V*)  "  Lih{t)  •  (7) 

Consequently,  if  the  assumption  of  linearity  over  a  time,  2T^,  can  be  made, 
the  standard  LI  correction  is  valid  using  sensors  of  voltage  and  current  at 
the  same  position. 
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Analogously,  to  the  extent  that  V±  and  I±  are  linear  from  t  -  T-  to  t  + 
T^,  equation  (5b)  reduces  to  1 

Ii  +  1(t)  =  ^(t)  -  CiVi(t)  ,  (8) 

where  is  the  capacitance  of  the  transmission  line. 


4.  IMPLEMENTATION  IN  VANDALS  PROGRAM  AT  AURORA 

VANDALS  is  a  program  implemented  on  the  Aurora1  PDP  1 1 /60  computer  to 
process  voltage  and  current  data  recorded  with  Tektronix  7912  transient  digit¬ 
izers  and  their  associated  software.  Although  VANDALS  performs  other  tasks 
associated  with  digitized  data-taking  (e.g.,  baseline  and  sweep-speed  adjust¬ 
ment?  time  alignment;  integration  of  traces  from  capacitive  voltage  dividers 
and  B-dot  loops;  and  calculation  of  wave  impedance,  power,  energy,  and 
IV** 2. 8),  its  main  function  is  to  implement  the  translation  algorithm  dis¬ 
cussed  above. 

The  relevant  portion  of  the  VANDALS  subroutine  that  performs  the  transla¬ 
tion  appears  in  appendix  A.  The  transmission-line  circuit  in  VANDALS  is  a 
string  of  up  to  30  transmission  lines  of  arbitrary  impedances,  Z(J),  and 
transit  times,  T(J).  At  Aurora  this  circuit  (fig.  2)  describes  the  region  of 
the  machine  from  the  high-voltage  end  of  the  vacuum  insulator  stack  to  the 
field-emission  load.  The  transit  times  of  the  circuit  elements  were  chosen  so 
that  the  various  voltage  and  current  monitors  and  the  load  are  positioned  at 
junctions  in  the  circuit. 

The  voltage  and  current  traces  to  be  translated  are  contained  in  the 
arrays  VI  and  II,  each  of  which  is  dimensioned  for  768  points.  These  data 
arrays  are  associated  with  any  one  (N)  of  the  31  positions  between  (or  at  the 
end  of)  transmission  lines  in  the  circuit.  Things  are  set  up  so  that  position 
N  is  between  lines  N  and  N  +  1.  Position  30  is  between  line  30  and  the  non¬ 
existent  line  31 ;  position  0  is  between  line  1  and  the  nonexistent  line  0. 
Each  pass  through  the  major  DO  loop  (ends  at  statement  200)  replaces  the  data 
in  the  arrays  VI  and  II  with  data  appropriate  to  either  position  N  +  1  or 
position  N  -  1.  The  formalism  of  the  algorithm  will  permit  translations  in 
either  direction  in  the  circuit  provided  an  adjustment  is  made  in  the  sign  of 
current  terms  in  the  statements  corresponding  to  equations  (5a)  and  (5b) 
above.  This  must  be  done  because  the  definition  of  positive  current  implies  a 
directionality  in  the  circuit.  The  factor  SIGN  which  takes  the  value  +1  or 
-1,  as  appropriate,  takes  care  of  this. 

As  explained  in  the  theory  section  above,  to  calculate  a  point  for  time  t 
in  a  translation  across  a  line,  J,  of  length  T(J),  data  must  exist  in  VI  and 
II  corresponding  to  both  t  +  T(J)  and  t  -  T(J).  This  means  that,  if  in  VI  and 
II  the  data  cover  the  time  span  from  t1  to  t2,  the  translated  data  (which  are 
temporarily  put  in  the  scratch  arrays  V2  and  12)  will  cover  the  time  span  from 

1B.  Bernstein  and  I.  D.  Smith,  Aurora,  and  Electron  Accelerator ,  IEEE  Trans . 
Nucl.  Sci.,  20_  (June  1973),  294. 
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Figure  2.  The  Aurora  sensor/load  data-correlation  circuit. 


t.|  +  T(J)  to  t2  -  T(J).  Therefore,  in  successive  translations  the  time  span 
covered  by  the  data  will  shrink.  The  program  keeps  track  of  this  with  the 
quantities  K1  =  NMN  and  K.2  =  NMX,  which  are  the  respective  numbers  of  the 
minimum  and  maximum  bins  of  relevant  data. 

Since  the  program  permits  transit  times  of  arbitrary  values,  an  interpola¬ 
tion  is  done  between  points  in  VI  and  II  to  get  quantities  V1(t±T(J))  and 
I 1 ( t±T ( J ) )  for  use  in  the  translation  formulas.  A  graphical  illustration  of 
this  interpolation  procedure  is  commented  into  the  listing  (app  A)  in  which 
all  parameters  of  relevance  are  identified.  For  our  purposes  a  linear  inter¬ 
polation  has  been  found  to  be  quite  adequate;  but,  as  a  result  of  this  inter¬ 
polation,  the  data  are  smoothed  in  translation.  This  smoothing  can  be  seen  in 
the  sample  Aurora  data  of  figure  3.  Also  shown  in  figure  3  is  the  shortening 
of  the  traces  due  to  the  algorithm.  This  shortening  is  at  both  ends,  although 
one  could  reason  that  it  need  not  be  for  times  before  the  pulse,  if  we  assume 
we  know  a  priori  what  the  current  and  voltage  are  before  the  time  frame  of  the 
recorded  data.  As  the  original  current  trace  can  be  seen  here  to  end  slightly 
before  the  voltage  trace,  in  this  example  it  is  the  extent  of  the  original 
current  data  that  governs  the  extent  of  data  in  subsequent  calculations. 


5.  TEST  OF  THE  TRANSLATION  ALGORITHM 

To  test  the  coding  of  the  basic  translation  algorithm,  a  simple  sine— wave 
example  is  presented,  with  a  very  simple  circuit. 

Consider  a  50-fl  transmission  line  that  is  100  ns  long  and  is  represented 
by  a  chain  of  twenty -five  50-ft  transmission  lines,  each  4  ns  long.  At  one  end 
of  this  line  some  unspecified  circuitry  is  attached.  This  circuitry  results 
in  a  voltage  of  50-V  amplitude  that  is  exactly  in  phase  with  the  current  of 
1  A  at  this  particular  end  of  the  transmission  line.  The  voltage  and  current 
are  sinusoids  with  a  250-ns  period. 

At  the  remote  end  we  expect  identical  waveforms  except  that  they  are 
delayed  by  100  ns.  The  load  responsible  for  V(t)  and  I(t)  must  be  a  resist¬ 
ance  of  50  (or  a  further  length  of  50-ft  cable  terminated  in  50  ft ) . 

Figure  4(a)  shows  the  given  voltage  and  current.  Figure  4(b)  shows  the 
voltage  and  current  calculated  for  the  other  end  of  the  1 00-ns  line.  Note 
expected,  the  traces  are  identical  in  shape  to  the  given  data  but 
delayed  by  100  ns.  Also  note  that  the  traces  are  shorter  by  about  120  ns  at 
each  end.  Theoretically,  the  data  should  be  shortened  by  one  transit  time  or 
100  ns  at  each  end.  The  extra  20  ns  is  due  to  a  rounding  effect:  on  each 
pass  through  one  of  the  25  transmission  lines,  the  data  are  shortened  by  an 
integral  number  of  data  points.  When  the  transit  time  is  4  ns  and  the  512 
data  points  are  0.98  ns  apart,  4.1  data  points  should  be  eliminated  per  trans¬ 
mission  line  at  each  end  of  the  trace.  To  round  4.1  down  to  4  data  points 
would  leave  an  irrelevant  point  at  each  end  of  each  trace.  Therefore,  5  data 
points  are  eliminated  instead;  (5  points)  x  (0.98  ns/data  point)  x  (25  lines) 
accounts  for  122  ns  lost  at  each  end  of  each  trace.  This  effect  has  proven  to 
be  no  problem  with  real  data  (fig.  3). 
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(a)  sensor  position 


Figure  3.  Sample  Aurora  data:  (a)  data  provided  by  the  sensors,  (b)  data 
calculated  for  the  load  position  with  the  use  of  the  data  of  3(a),  and 
(c)  data  calculated  for  the  sensor  position  based  on  the  calculated  load- 
position  data  of  3(b).  Ideally  this  would  reproduce  exactly  the  data  in 
3(a),  but  is  smoothed  due  to  interpolation. 
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6.  CONCLUDING  REMARKS 


Legitimately,  one  can  only  make  use  of  the  equations  (5a)  and  (5b)  if  V 
and  I  are  known  functions  of  time  at  some  common  position.  However,  if  expe¬ 
rience  indicates  that  either  V  or  I  does  not  vary  too  much  in  position,  one 
could  rationalize  applying  this  method  to  two  traces  from  slightly  different 
positions  along  a  transmission  line.  The  present  method  of  analysis,  which 
assumes  that  the  two  traces  originate  at  the  same  position  along  a  transmis¬ 
sion  line,  results  in  a  closed  form  (equations  (5a)  and  (5b))  for  translating 
the  traces.  One  could  envision  converging  infinite  series  replacements  for 
equations  (5a)  and  (5b)  so  that  information  for  a  third  position  could  be 
unfolded  rigorously  from  monitors  at  two  different  positions.  An  approach  to 
this  problem  was  taken  by  R.  Richardson  et  al.  of  Maxwell  Laboratories,  Inc. 
in  reference  to  work  on  their  Blackjack  5  Machine.2 


Figure  4.  Sinusoidal  test  of  code: 
(a)  data  simulated  for  "sensor" 
position  and  (b)  data  calculated 
(using  VANDALS)  100  ns  downstream 
of  original  data  through  twenty- 
five  4-ns  transmission  lines.  The 
extra  time  missing  at  each  end  of 
the  calculated  data  results  from 
rounding  up  the  fractional  number 
of  bins  that  should  be  removed  from 
the  data  on  each  pass. 


(a)  original  data 


(b)  after  translation 


CR.  Richardson,  E.  Chu,  W.  Clark,  J.  Shannon,  and  M.  Wilkinson,  Calibration 
of  the  Blackjack  5  Pulse  Generator  Output  Power,  National  Bureau  of  Standards 
Special  Publication  628,  Measurement  of  Electrical  Quantities  in  Pulse  Power 
Systems,  R.  H.  McKnight  and  R.  E.  Hebner,  Jr .,  ed.(June  1982),  150. 
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APPENDIX  A 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


MPLUS+1 


MMINUS+1 


MMINUS 


FIRST  TBTUN- 


MPLUS 


-SECOND  TBTUIN 


TIMEX) 


MMINUS 
OR  — 
I MINUS 


•>X 


J  X< 


URL  US 
— OR 
I  PLUS 


1 


9 


•TMINUS- 


* 

3 


-+- 

4 


6 


■T  (  J  )  • 


7 


T(  J)- 


9 


•-+- 

10 

* 


:l  l 


-K*DT- 


•TPL  US  • 


THE  GRAPHIC  EXAMPLE  SHOWS  EITHER  A  VOLTAGE  OR  CURRENT  DATA  ARRAY 
THAT  IS  TO  CONTRIBUTE  TO  A  DATA  TRANSLATION ♦  THE  NUMBERS  IDENTIFY 
THE  BINS ♦  EACH  *  IS  A  DATA  POINT  IN  THAT  ARRAY ♦  THE  SYMBOL 
X  SHOWS  THE  LINEARLY  INTERPOLATED  VALUES  USED  IN  THE  ALGORITHM 
(STATEMENT  100  AND  THE  ONE  JUST  ABOVE  IT)  THAT  ACCOMPLISHES  THE 
TRANSLATION ♦ 


K 1  AND  K2  ARE 


THE  respective  scratch  equivalents  of  nmn  and 


NMX. 


300  K 1 "  NMN 
K2--  NMX 

THE  SCRATCH  ARRAYS  V2  AND  12  ARE  ZEROED* 

DO  99  K=  1 t  768 
V2(K>=  0. 

99  I 2 ( K ) ~  0* 
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c 

c 

c 

e 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


POSITIVE  CURRENT  IS  ASSUMED  TO  BE  IN  THE  DIRECT TON 
VALUES  FOR  THE  SUBSCRIPT  J  IN  Z(J> 


USED  AS  THE  SIGN  OF  THE  CURRENT  TERMS 


OF  INCREASING 

AND  T < J ) ♦  SIGN  IS  A  FACTOR  +  OR  ~l 
I N  T  H  E  T  R  A  N  S  L  A  T 1 0  N  A  L  G  0  R I T  H  M  * 

TO  USE  THE  ALGORITHM  FOR  TRANSLATIONS  IN  THE  DIRECTION  OF  INCREASING 

nrnrrrTni8  tlU ,T?oUSf  THE  ALGORITHM  FOR  TRANSLATIONS  IN  THE  OTHER 
l-'ihtL  I  lUNf  a J.l)N  lb  —  1  , 


SIGN=  1 , 

IF<N1  ♦ GT .  N2)  SIGN-  -1, 

IF  (N1  ,EQ.  N2)  GO  TO  1000 
NIJMBER=  ABS  <  N2-N1 ) 

THE  LOOP  OVER  CIRCUIT  ELEMENTS  BEGINS. 


DO  200  JJ=1 ?  NUMBER 


J  IS  THE  I.D. 
IMPLICIT  HERE 


NUMBER  OFT  H  E  E LEM E N  T  B E I N  G  CAL C U L  A  T E D  A  C  R  0 S S . 

1'S  THAT  TO  THE  RIGHT  OF  POSITION  J  IS  ELEMENT  J+.I.  . 


J=  JJ  +  N:l. 

I P ( N 1  , GT .  N2 )  J=N1+1~JJ 


DT  IS  THE  TIME  PER  BIN. 

THE  IDENTIFIERS  K1  AND  K2  OF  THE  MINIMUM  AND  MAXIMUM  BINS 


DATA  ARE  ADJUSTED  BELOW.  IF  T1  (T2)  IS  THE  TIME  AT  WHICH 
BEGINS(ENBS)  BEFORE  TRANSLATION  i  Tl  +  T  ( T2-T )  IS  THE  TIME 
THE  DATA  BEGINS  (ENDS)  IN  THE  NEW  TRANSLATED  DATA  ARRAY. 


OF  RELEVANT 
THE  DATA 
AT  WHICH 


K 1 ~  1  +  Kl+T  < J)/DT 
K2-  K2-T < J)/DT 
I FLAG  =  0 

IF  (K2  .GT.  Kl)  GO  TO  10 
WRITE  (5,  11) 

11  FORMAT  (  ■'  THE  SPAN  OF  THE  DATA  IS  TOO  RFSTR  ITTFD 
1  '  TRANSLATION.') 

I FLAG  =  1 
GO  TO  1000 


FOR  FURTHER ' > 


K.l  AND  K2  ARE  RESPECTIVELY  THE  LOWEST  AND  HIGHEST  BINS  OF  RFI  FVANT 
DATA.  THE  LOOP  OVER  BINS  IN  THE  CALCULATED  DATA  ARRAYS  BEGINS  HERE, 
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DO  100 

TMINUS- 

MMINUS- 

TBTWN- 

VMINUS= 

IMINUS:: 

TPLUS- 

MPLl.JS= 

TBTWN-- 

VPLUS= 

IPL.US= 


K~  K.l  y  K2 
=  KYDT-T ( J ) 

=  TMINUS/DT 
TMINUS  ~  MM IN US  *  DT 

:  VI  (MMINUS)-!TBTWN*(V1  (MMINUS+1  )-Vl  (MMIN(JS)  ) /DT 
K>kl  )  +TBTWN#  (  1 .1.  ( MMINUS+1 )  -  1 1  ( MMINI.JS)  )  /DT 

TPLUS/DT 
TPLUS-MPLUS+D T 

VI  ( MPLUS )  +TBTWN*  ( V:L  ( MPLUS  +  l ) ~V1  ( MPLUS )  ) /DT 
I 1 ( MPLUS ) +TBTWN* ( I 1 ( MPLUS+l ) -I 1 ( MPLUS ) ) /DT 
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no  non  non o 


APPENDIX  A 


C  THE  NEXT  TWO  STATEMENTS  DO  THE  ACTUAL  TRANSLATION « 

V2  <  K ) =  , 5*(VMINUS+SIGN*Z( J > #IMINUS+VPLUS-SIGN*Z < J)# I PLUS) 

100  I2(K)~  .5#<SIGN*VMINUS+Z< J>*IMINUS-SIGN#VPLUS+Z< J>*IPLUS)/Z< J> 

THE  TRANSLATED  DATA  IN  THE  SCRATCH  ARRAYS  ARE  TRANSFERRED  TO  THE 
ORIGINAL  ARRAYS  AND  THE  SCRATCH  ARRAYS  ARE  ZEROED* 

DO  200  K-l r  768 
01 (K>=02<K) 

02  <  K ) =0 . 

I1(K)-I2(K) 

200  12  <  K ) “0 ♦ 

THE  LOWEST  AND  HIGHEST  RELEVANT  DATA  BINS  (NMN  AND  NMX)  ARE  SHIFTED 
SINCE  DATA  IS  LOST  AT  EACH  END  OF  THE  TRACE  BY  AN  AMOUNT  CORRESPONDING 
TO  THE  TRANSIT  TIME  OF  THE  TRANSLATION. 

NMX=  K2 
NMN=K1 
1000  RETURN 
END 
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